Decoherence by a spin thermal bath: Role of the spin-spin interactions 
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We study the decoherence of two coupled spins that interact with a spin-bath environment. It is 
shown that the connectivity and the coupling strength between the spins in the environment are of 
crucial importance for the decoherence of the central system. Changing the connectivity or coupling 
strenghts changes the decoherence of the central system from Gaussian to exponential decay law. 
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It is commonly accepted that decoherence by nuclear 
spins is the main obstacle for realization of quantum com- 
putations in magnetic systems; see, e.g., discussions of 
specific silicon [l[ and carbon Q based quantum com- 
puters. Therefore, understanding the decoherence in 
quantum spin systems is a subject of numerous works 
(for reviews, see Refs 0, The issue seems to be 

very complicated and despite many efforts, even some 
basic questions about character of the decoherence pro- 
cess are unsolved yet. Due to the interactions with and 
between the spin of the bath, an analytical treatment 
can be carried out in very exceptional cases, even if the 
central systems contains one spin only. Recent work 
suggests that the internal dynamics of the environment 
can be crucial to the decoherence of the central sys- 
tem @, H, 0, 0, I, EH EH, d El EI, 
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Letter, we present results of extensive simulation work 
of a two-spin system interacting with a spin-bath envi- 
ronment and show that the decoherence of the two-psin 
system can exhibit different behavior, depending on the 
characteristics of the coupling with the environment and 
of the internal dynamics of the latter. We also provide a 
simple physical picture to understand this behavior. 

In general, the behavior of an open quantum system 
crucially depends on the ratio of typical energy differ- 
ences of the central system 5E C and the energy E ce which 
characterizes the interaction of the central system with 
the environment. The case 6E C <C E ce has been studied 
extensively in relation to the "Schrodinger cat" problem 
and the physics is quite clear E3|: As a result of 
time evolution, the central system passes to one of the 
"pointer states" [l^l which, in this case, are the eigen- 
states of the interaction Hamiltonian H ce . The oppo- 
site case, 6E C ^> E ce is less well understood. There is 
a conjecture that in this case the pointer states should 
be eigenstates of the Hamiltonian H c of the central sys- 
tem but this has been proven for a very simple model 
only [20| . On the other hand, this case is of primary in- 
terest if, say, the central system consists of electron spins 
whereas the environment are nuclear spins, for instance 
if one considers the possibility of quantum computation 
using molecular magnets 2lJ, |22j . 



We consider a generic quantum spin model described 
by the Hamiltonian H = H c + H ce + H e where H c = 
— JSi ■ S2 is the Hamiltonian of the central system and 
the Hamiltonians of the environment and the interaction 
of the central system with the environment are given by 
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respectively. The exchange integrals J and de- 
termine the strength of the interaction between spins 
S„ = (S%,Sv,Sn) of the central system, and the spins 
I„ = (I*, 1^,1*) in the environment, respectively. The 
exchange integrals A.\ a ) control the interaction of the cen- 
tral system with its environment. In Eq. (frj, the sum 
over a runs over the x, y and z components of spin- 
1/2 operators S and I. In the sequel, we will use the 
term "Heisenberg-like" H e (H ce ) to indicate that each 
Q 



(A^p) is a uniform random number in the range 
[— f2|,f2] ([—A, AD, and A being free parameters. In 
earlier work lfl 12 1 , we found that a Heisenberg-like H e 
can induce close to perfect decoherence of the central 
system and therefore, we will focus on this case only. 

The bath is further characterized by the number of en- 
vironment spins K with which a spin in the environment 
interacts. If K = 0, each spin in the environment in- 
teracts with the central system only. K = 2, K = 4 or 
K = 6 correspond to environments in which the spins are 
placed on a ring, square or triangular lattice, respectively 
and interact with nearest-neighbors only. If K — N — 1 , 
each spin in the environment interacts with all the other 
spins in the environment and, to give this case a name, 
we will refer to this case as "spin glass" . 

If the Hamiltonian of the central system H c is a per- 
turbation relative to the interaction Hamiltonian H cei 
the pointer states are eigenstates of H ce [lj|. In the 
opposite case, that is the regime |A| -C |J| that we 
explore in this Letter, the pointer states are conjec- 
tured to be eigenstates of H c The latter are given 
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by |1) = |T0 = |TT), |2> = \S) = (|U) - UT))/V2, 
|3> = |T ) = (|TI> + UT))/V2, and |4) ee |T_r) = |U>, 
satisfying H C \S) = (3J/4)|S) and H c \T t ) = (-J/4)|T;) 
for i = -1,0, 1. 

The simulation procedure is as follows. We generate 
a random superposition \<j>) of all the basis states of the 
environment. This state corresponds to the equilibrium 
density matrix of the environment at infinite tempera- 
ture. The spin-up - spin-down state (|TI)) is taken as 
the initial state of the central system. Thus, the initial 
state of the whole system reads \9(t — 0))) = |T|) and 
is a product state of the state of the central system and 
the random state of the environment which, in general is 
a (very complicated) linear combination of the 2 N basis 
states of the environment. In our simulations we take 
N = 16 which, from earlier work 1(| 17 1 , is sufficiently 
large for the environment to behave as a "large" system. 
For a given, fixed set of model parameters, the time evo- 
lution of the whole system is obtained by solving the 
time-dependent Schrodinger equation for the many-body 
wave function describing the central system plus 

the environment [23| . It conserves the energy of the whole 
system to machine precision. We monitor the effects of 
the decoherence by computing the the matrix elements 
of the reduced density matrix p (t) of the central system. 
As explained earlier, in the regime of interest |A| <§; \ J\, 
the pointer states are expected to be the eigenstates of 
the central systems. Hence we compute the matrix ele- 
ments of the density matrix in the basis of eigenvectors 
of the central system. We also compute the time depen- 
dence of quadratic entropy S c (t) = 1 — Trp 2 (t) and the 
Loschmidt echo L (t) = Tr (p (t) p (t)) [H, where p (t) 
is the density matrix for H ce = 0. 

In Fig. [1] and Fig.[2j we show the time evolution of the 
elements of the reduced density matrix p (i) for different 
connectivity K and Q, for the case that H ce is an isotropic 
Heisenberg model. 



If |A| ^> KQ, in agreement with earlier work [25l. |26|. 
we find that in the absence of interactions between the 
environment spins (KQ = 0) and after the initial de- 
cay, the central system exhibits long-time oscillations (see 
Fig. QIa)(left)). In this case and in the limit of a large 
environment, we have [2f| 



Re p 23 (t) 
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where b = NA 2 /4, c = b/2 and uo = J - A. Equation © 
clearly shows the two-step process, that is, after the ini- 
tial Gaussian decay of the amplitude of the oscillations, 
the oscillations revive and their amplitude levels of [26| . 
Due to conservation laws, this behavior does not change if 
the environment consists of an isotropic Heisenberg sys- 
tem = Q for all a, i and j), independent of K. 
If, as in Ref. [25|, we take A 
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FIG. 1: (Color online) The time evolution of the real part of 
the off-diagonal element pn (left panel) and the diagonal ele- 
ments pn, . . . , /944 (right panel) of the reduced density matrix 
of a central system (with J = —5), coupled via an isotropic 
Heisenberg interaction H ce (A = —0.075 ) to a Heisenberg- 
like environment H e (Q — 0.15) with different connectivity: 
(a) K = 0; (b) K = 2; (c) K = 4; (d) K = 6; (e) K = N - 1. 



long-living oscillations is no longer constant but decays 
very slowly [25|. 

If |A| w Kil, the presence of Heisenberg-like inter- 
actions between the spins of the environment has little 
effect on the initial Gaussian decay of the central sys- 
tem, but it leads to a reduction and to a decay of the 
amplitude of the long-living oscillations. The larger K 
(see Fig. QIb-e)(left)) or fi (see Fig. H>,c)), the faster 
the decay is. Note that for the sake of clarity, we have 
suppressed the fast oscillations by plotting instead of the 
real part, the absolute value of the matrix elements. 

If |A| <C KQ, keeping K fixed and increasing Q 
smoothly changes the initial decay from Gaussian (fast) 
to exponential (slow) , and the long-living oscillations are 
completely suppressed (see Fig.[UJb,d)). For large Q, the 
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FIG. 2: (Color online) The time evolution of the off-diagonal 
element P23 of the reduced density matrix of a central system 
(with J = —5), interacting with a Heisenberg-like environ- 
ment H s via an isotropic Heisenberg Hamiltonian H ce (with 
A = —0.075 ) for the same geometric structures in the envi- 
ronment: (a,b) K = 2 and (c,d) K — N — 1. The number 
next to each curve is the corresponding value of Q. 



(b) K=2 




100 120 140 160 180 200 
t 



FIG. 3: (Color online) Same as Fig. [2] except that H ce is 
Heisenberg-like and A = 0.15. 



simulation data fits very well to 

\P23 (*)| 



1 -A K (n)t 



(3) 



with A K (ft) « flA K , A 2 = 9.13 and A N - X = 26.73. 

Physically, the observed behavior can be understood as 
follows. If |A| ss K£l, a bath spin is affected by roughly 
the same amount by the motion of both the other bath 
spins and by the two central spins. Therefore, each bath 
spin has enough freedom to follow the original dynamics, 
much as if there were no coupling between bath spins. 
This explains why the initial Gaussian decay is insensitive 
to the values of K or ft. After the initial decay, the 



whole system is expected to reach an stationary state, but 
because of the presence of Heisenberg-like interactions 
between the bath spins, a new stationary state of the bath 
is established, suppressing the long-living oscillations. 

For increasing K, the distance between two bath spins, 
defined as the minimum number of bonds connecting the 
two spins, becomes smaller. For instance, for K = 2, 
this distance is (N — 2) /2, and for K = N — 1, it is zero. 
Therefore, for fixed fl and increasing K the fluctuations 
in the spin bath can propagate faster and the evolution to 
the stationary state will be faster. Similarly, for fixed K, 
increasing the coupling strength between the bath spins 
will speed up the dynamics of the bath, that is, the larger 

the faster will be the evolution to the stationary state. 

In the opposite case |A| <C K£l, H ce is a small per- 
turbation relative to H e and the coupling between bath 
spins is the dominant factor in determining the dynamics 
of the bath spins. Therefore, by increasing K or O, the 
bath spin will have less freedom to follow the dynamics 
induced by the coupling to the two central spins, the in- 
fluence of the bath on the central system will decrease, 
and the (exponential) decay will become slower. 

According to the general picture of decoherence [lij ]. 
for an environment with nontrivial internal dynamics 
that initially is in a random superposition of all its eigen- 
states, we expect that the central system will evolve to 
a stable mixture of its eigenstates. In other words, the 
decoherence will cause all the off-diagonal elements of 
the reduced density matrix to vanish with time. In the 
case of an isotropic Heisenberg coupling between the cen- 
tral system and the environment, H c commutes with the 
Hamiltonian H, hence the energy of the central system is 
a conserved quantity. Therefore, the weight of the singlet 
IS*) in the mixed state should be a constant (1/2), and 
the weights of the degenerate eigenstates \T ), |T-i) and 
|7i) are expected to become the same (1/6). As shown 
in Fig. [Qb-e) (right), our simulations confirm that this 
picture is correct in all respects. 

It is important to note that although in the foregoing 
discussion we have compared KQ, to |A|, this does not 
imply that KQ can be used to fully characterize the de- 
coherence process. In order to clarify the role of K and 
f2, we change the coupling between the central system 
and the bath from Heisenberg to Heisenberg-like. From 
a comparison of the data in Fig. [2] and Fig. [3l it is clear 
that the roles of K and are the same in both cases, no 
matter whether the central-bath coupling is isotropic or 
anisotropic. However, there are some differences in the 
decoherence process. 

If |A| ^> Kfl, in the presence of anisotropic interac- 
tions between the central system and the environment 
spins, even in the absence of interactions between the 
bath spins, the second step of the oscillations decays and 
finally disappear as K increases. This is because the 
anisotropic interactions break the rotational symmetry 
of the coupling between central system and environment 



4 




herence [19|. Furthermore, we have demonstrated that, 
in the case that the environment is a spin system, the 
details of this spin system are important for the decoher- 
ence of the central system. In particular, we have shown 
that changing the internal dynamics of the environment 
(geometric structure or exchange couplings) may change 
the decoherence of the central spin system from Gaussian 
to exponential decay. 
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FIG. 4: (Color online) The time evolution of the the entropy 
S c (t) and Loschmidt echo L (t) of a central system (with 
J = —5), interacting with a Heisenberg-like environment H e 
(with different Q) via a Heisenberg (a,b, A = —0.075) or 
Heisenberg-like (c,d, A = 0.15) Hamiltonian H C£ for the case 
K = 2. The number next to each curve is the corresponding 
value of Q. 

which is required for the long-living oscillations to per- 
sist. 

If |A| <C Kfl, \p 23 (t)| can still be described by Eq. ©, 
but now Ak (Q) is no longer a linear function of fi. For 
an anisotropic H ce , the energy of the central system is 
no longer a conserved quantity. Therefore there will be 
energy transfer between the central system and the envi- 
ronment and the weight of each pointer state (eigenstate) 
in the final stable mixture need not be the same for all 
K or n. 

For a change, we illustrate this point by considering 
the quadratic entropy S c (t) and Loschmidt echo L (t) . 
We expect that these quantities will also dependent of 
the symmetry of the coupling between central system 
and the spin bath. In Fig. [4J we present results for 
large and K = 2, confirming this expectation. For 
isotropic (Heisenberg) H ce and perfect decoherence (zero 
off-diagonal terms in the reduced density matrix) we ex- 
pect that m&x t S c (t) = 1- [(l/2) 2 + 3x (1/6) 2 ] = 2/3, in 
concert with the data of Fig. ID^a)). For Heisenberg-like 
H ce , max ( S c (t) will depend on the coupling strengths 
and as shown in Fig. [He), we find that max t 5' c (i) = 
1 — 4 x (1/4) = 3/4, corresponding to the case that all 
the diagonal elements in the reduced density matrix are 
the same (1/4) and all other elements are zero. 

In conclusion, we have shown how a pure quantum 
state of the central spin system evolves into a mixed state, 
and that if the interaction between the central system 
and environment is much smaller than the coupling be- 
tween the spins in the central system, the pointer states 
are the eigenstates of the central system. Both these ob- 
servations are in concert with the general picture of deco- 
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